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Abstract 

The ratio of the off-diagonal hyperfine amplitude to the tensor transition 

polarizability (M^f/ '/?) for the 6S— >7S transition in cesium has been measured. 

The value of f3 = 27.024(43) exp t(67)theoryao * s then obtained using an accurate 

semi-empirical value of Mm. This is combined with a previous measurement 

of parity nonconservation in atomic cesium and previous atomic structure 

calculations to determine the value of the weak charge. The uncertainties in 

the atomic structure calculations are updated (and reduced) in light of new 

experimental tests. The result Qw = — 72.06(28) CX pt(34)theory differs from the 

prediction of the standard model of elementary particle physics by 2.5a. 
PACS number(s):32.80.Ys, 11.30.Er, 12.15.Ji, 32.10.Dk 
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I. INTRODUCTION 



Electroweak experiments have now reached high precision in testing the standard model 
and in searching for new physics beyond it . These experiments include measurements of 
parity nonconservation (PNC) in atoms as first proposed in Ref. 0. Atomic PNC measure- 
ments are uniquely sensitive to a variety of new physics, such as the existence of additional 
Z bosons, because of the different energy scale and because they probe a different set of 
model-independent electron-quark coupling constants than those measured by high-energy 
experiments [[J. The most precise atomic PNC experiment M examines the mixing of S 
and P states in atomic cesium. Specifically, it compares the mixing due to the PNC neutral 
weak current interaction to the S-P mixing caused by an applied electric field ( "Stark mix- 
ing"). In previous work [|J, this measurement was combined with theoretical calculations 
of the structure of the cesium atom to obtain the weak charge Qw, which characterizes the 
strength of the neutral weak interaction and can be compared to the value predicted by 
the standard model. The atomic structure calculations were used to obtain two pieces of 
information: the amount of Stark mixing and the relevant PNC electronic matrix elements. 
The 1.2% uncertainty in the determination of Qw was dominated by the uncertainties in 
those two calculated quantities. In this paper we report a reduced uncertainty in Q w that is 
obtained by 1) measuring the Stark mixing, and 2) incorporating new experimental data into 
the evaluation of the uncertainty in the calculation of the PNC matrix elements. These new 
data indicate that the calculations are more accurate than was indicated by the less precise 
(and in some cases incorrect) data available at the time the calculations were published. 



II. THEORY 

The QS ground state and 7S excited state of atomic cesium both have two hyperfine 
levels: F = 3 and F = 4. In the presence of a dc electric field E, a magnetic field, and 
a standing- wave laser field with propagation vector k and polarization e, the AF = ±1 
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6S — > 7S amplitudes, used in both Ref. [|J and the present work, are given by 

A 6S ^ 7S = [i(3{E x e ) + Ml(k x e) + El PNC e] ■ (F'm' F \ a \ Fm F ), (1) 

where Ml = M ± Mhf5_FF'±i is the magnetic dipole amplitude (M is from relativistic and 
spin-orbit effects, is from the off-diagonal hyperfine interaction), and a is the Pauli 
spin matrix. The tensor transition polarizability || (3 characterizes the size of the Stark 
mixing-induced electric dipole amplitude, and E'Ipnc is the PNC matrix element given by 



Q 



E1 PNC = (75 \D\ 6S) = ^ fc PNC . (2) 



Here, | nS) is an | nS) state into which the PNC Hamiltonian has mixed a small amount of 
| nP) states, D is the electric dipole operator, N is the number of neutrons, and /cpnc is the 
calculation of the sum of relevant matrix elements between S and P states given by 
N ((7S\D\ nP){nP \ H PNC I 6S) , (7S \ H PNC \ nP)(nP \ D \ QS}\ 

Since ifpNc = Gf^QwPn{^) I \/&, each of the terms in Eq. (|3|) is the product of a dipole 
matrix element times a 7 5 matrix element evaluated in the nucleus. Ninety-eight percent of 
the sum comes from the 6Pi/ 2 and 7P±/2 states [|]. 

In Ref. [|J, Im(ElpNc)//? is measured. The value Qw is obtained by multiplying this 
ratio by fiN/kp^c- This paper concerns the improved determination of (3 and fcpNc, and 
thus Q w . 

To determine (3, we measure Mj,f//3 and take advantage of the fact that Mm can be 
accurately determined semi-empirically [0. The amplitude M hf is due to the hyperfine 
interaction and thus can be expressed in terms of well-measured hyperfine splittings. In 
this experiment we observe the QS — > 7S rate driven with a standing-wave laser beam with 
polarization e = ez and a field geometry (E along x) such that the transition rate is 

| A 6S ^ 7S | 2 = (3 2 E 2 e 2 + (M ± M hf <Wi) 2 e 2 , (4) 

where small interference terms have been omitted. The /3-PNC and Ml-PNC interference 
terms are negligible and the (3-Ml interference terms cancel almost identically (< 10~ 6 ) 



because of their k dependence and the standing- wave geometry of the experiment. We 
determine M^/(3 by measuring the total rate on the two AF = ±1 hyperfme transitions 
with large E, where | A 6S ^ 7S | 2 ~ f3 2 E 2 , and with E = 0, where | A 6 s—7S | 2 ~ {M±M h{ 5FF>±i) 2 - 
We combine the ratios of the high and low E rates on both transitions to determine M^/ j3. 

A complication arises because the locations of the antinodes of the oscillating electric 
(e ac ) and magnetic (b ac ) fields are separated by A/4 in the standing wave. Because of this 
separation, photoionization (which is driven by e ac ) is larger for 7S atoms excited by e ac (El 
atoms) than it is for 7S atoms excited by b ac (Ml atoms). The result is that the detection 
efficiency for El excitations is slightly smaller (~ 1% for typical intensities) than for Ml 
excitations. This difference gives a potential systematic error that is intensity dependent. 
The ratios of the signals, measured at a laser intensity /, for the AF = +1 and AF = — 1 
transitions, respectively, are then 

ir 4 = (^Z^) 2 (l+r,/),and (5a) 

^M^flUrt, (5b) 
where rj is a parameter that describes the difference in photoionization fraction. 



III. EXPERIMENT 

The apparatus used in the present experiment is very similar to that in Refs. [f|||. A 
collimated beam of cesium is optically pumped into either the F = 3 or F = 4 hyperfine 
level of the QS1/2 ground state. The beam of atoms then travels roughly along the z axis 
into a region with mutually orthogonal dc electric (along x) and magnetic (along z) fields 
and intersects a 540 nm standing-wave laser field (along y) at right angles. The laser field 
is produced by a tunable dye laser that is frequency-locked to a finesse ~ 10 5 Fabry-Perot 
etalon. The etalon is, in turn, locked to a stable reference cavity. The light going to 
the reference cavity is double-passed through an acousto-optic modulator (AOM), so the 
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frequency of the laser light interacting with the atomic beam is z/i aser = Reference — 2i/aom- 
Thus, we can change the frequency of the dye laser in a very controlled manner by changing 
the frequency of the AOM. The dye laser drives the 65* — > 75 transition. Approximately 
half of the atoms excited to the 75 state relax to the previously depleted hyperfme ground 
state (F = 3 or F = 4). Further downstream, the atoms in the repopulated hyperfme 
level scatter photons from a diode laser probe beam tuned to an appropriate 6S1/2-6P3/2 
cycling transition. We collect the scattered photons on a large-area photodiode, and its 
photocurrent is proportional to the number of atoms making the 6S— >7S transition. 

To measure the ratio R^ A (or Rj^ 3 ), we scan the laser over the 65 — ► 75 AF = +1 
(or AF = —1) transition in 0.3 MHz steps. After each step we integrate the photocurrent 
for 16.67 ms and store that data point on disk. We alternate between scans with E = 
707.63(68) V/cm and E = V/cm. 

There is a 540 nm-laser-frequency-independent background signal from atoms in the 
wrong hyperfme state that is ~ 100 times larger than the desired Ml signal for E — V/cm. 
We measure this background before and after each data point by detuning the laser ~ 50 MHz 
from line center and measuring the photocurrent. These background points are measured 
alternately above and below the line center to cancel any linear frequency dependence of the 
background. We subtract the average background from the data points to leave only the 
contribution from atoms making the 6S— >7S transition. The sum of all the data points (the 
area under the spectral line) is proportional to the total transition rate. 

We looked for but did not observe any frequency dependence to the background. Also, 
all likely mechanisms, such as molecular transitions or light scattering off the mirrors, should 
have very broad spectral features, and hence will be eliminated by the background subtrac- 
tion. The uncertainty in our results due to possible frequency dependent backgrounds is less 
than 0.05%. 

Sample background-subtracted scans are shown in Fig. [I]. The two line shapes are asym- 
metric and slightly offset from one another because of their differing sensitivity to ac Stark 
shifts as discussed in Ref. ||. The different line shapes do not affect our measurement of 
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the total transition rate because the atoms' total transition amplitude is unchanged, even 
though the resonant frequency of each atom is shifted according to the local e ac field. There- 
fore, by integrating the areas under the entire broadened lines we can determine the desired 
relative ratios R^ 4 and Rj^ 3 . 



IV. RESULTS 

The detection efficiency and signal-to-noise ratio are significantly higher for i? 3 ^ 4 ; we 
measure that ratio at five different intensities from 0.6 kW to 2.8 kW and determine r\ 
to 1.5 parts in 10 3 using a least squares fit. We find the ratios Rq~* a = 2.4636(8) x 10~ 3 
and R A ~* 3 = 1.1357(6) x 10~ 3 where the two uncertainties have a common contribution 
from the extrapolation to zero intensity. Combining these results using Eq. (|5|) we find 
M M /(3 = -5.6195(91) V/cm. 0. From Ref. we take M hf = -151.86(38) (V/cm)a 3 , 



which is based on measured hyperfine splittings with a 0.3 ± 0.3% theory correction due to 
many body effects. This gives 

(3 = 27.024(43) cxpt (67) thcory a 3 . (6) 

This value is in excellent agreement with the semi-empirical values (3 = 27.17(35)aj] |7J and 
(3 = 27.15(13)a 3 Q, and the calculated value (3 = 27.00a 3 Q. 



Using our measured values for (3 and Im(i£lpNc)/A and the calculated value of fcpNc, 
we can now extract Qw The key issue is the uncertainty in the value of /cpnc- The au- 
thors of Refs. ||[r^-[14| discuss this issue at considerable length. Here we only summarize 
the conclusion of both groups that the most reliable measure is to use the same ab initio 
calculations of the electronic structure that are used to find fcpNc to calculate dipole ma- 
trix elements and hyperfine splittings for the 65*1/2, 7Si/2, 6P1/2, and 7Pi/ 2 states. The 
differences between these calculated values and the experimental determinations provide a 
reliable quantitative indication of the uncertainties in the calculations of /cpnc- The authors 
considered how well these errors in the hyperfine splittings and dipole matrix elements re- 
flect errors in /cpnc by rescaling their calculations in a variety of ways and comparing the 



relative sensitivities of the different quantities. They found that /cpnc has comparable or 



smaller sensitivity than the other quantities [15|. From comparing calculated and measured 



quantities, both groups arrived at uncertainties of about 1% for their value of fcpNc- Since 
the time that Refs. P,[1"2}-|T4]| were published, there have been a number of new and more 



precise measurements of the quantities of interest. In all cases, the new measurements show 
better agreement with the calculations than earlier measurements and also show that the 
largest previous disagreements were likely due to experimental errors. 

In Table [j] we have collected the results of the most precise measurements of relevant 
quantities in cesium. We list the quantities measured, the primary aspect of the electronic 
wave functions that is being tested in each comparison, and the difference between theory 
and experiment. Particularly notable are the top three lines of the table, which show that the 
agreement has dramatically improved from the 1-2% disagreements of the older experiments. 
In addition to the data in this table, there have been new experiments that revealed errors 
in earlier lifetime measurements in sodium and lithium. These new data eliminate what had 
appeared to be troubling 1% errors in equivalent calculations for those atoms. 

The standard deviation of the fractional differences between theory and experiment in 
Table III is 4.0 x 10 -3 . We believe this to be the most valid number to use to represent 



the 68% confidence level for fcpNc- Using the average of fcpNc = 0.905 x 10 n ieao fl2| and 



/cpnc = 0.908 x 10 n iea [[L3[] this gives a value of /c PNC = 0.9065(36) x 10 u iea 



When combined with our new value for f3 and the experimental PNC measurement this 
gives 

Q w = -72.06(28) cxpt (34) theory . (7) 

The standard model value including radiative corrections is Qw = —73.20(13) [HJ. Adding 
the uncertainties in quadrature, these values differ by 2.5 a. 

Assuming that this difference is not due to an experimental error or a statistical fluctua- 
tion, it suggests several possibilities. The first possibility is that the calculated value of the 
7 5 matrix element is in error by the requisite 1.58%. In light of Table I, such an error would 



require a wave function with a somewhat peculiar and insidious shape. Although none of 
the measured quantities depends on the shape of the wave function in a manner identical to 
that of 7 5 , the different comparisons in Table I do probe the value of the wave function in 
all regions: short, intermediate, and long distances. The largest single difference of the 16 
comparisons is only 0.79%, and the standard deviation is only 0.40%. The second possibil- 
ity is that there are contributions or corrections to atomic PNC within the standard model 
that have been overlooked. We see no justification for either of these two possibilities, but 
they clearly need to be explored further. The first offers a formidable but not overwhelming 
challenge to both theoretical and experimental atomic physicists. 

The final possibility is that this discrepancy is indicating the presence of some new 
physics not contained in the standard model. Physics that would be characterized by the 



S parameter |T7 is not a likely candidate because the size of the contribution needed [S = 
— 1.4(6)] would be in conflict with other data [l|]. However, there are other types of new 
physics, such as an additional Z boson, that would be consistent all other current data. 
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TABLE I. Fractional differences (xlO 3 ) between measured and calculated values of quantities 
relevant for testing PNC calculations in atomic cesium. We only list the most precise experiments. 
The second column lists the most relevant aspects of the wavefunctions that are being tested. 
(l/r 3 ) n p is the average of 1/r 3 over the wavefunction of the electronic state nP. Where the 
experiment has improved or changed significantly since the publication of Ref. (T^j, the difference 
from the old experiment is listed in brackets. 



Quantity 


Calculation 


Difference (xlO 3 ) 




measured 


tested 


Dzuba, et a/. a,b Blundell et al. c 


&Expt 



6S^7S dc Stark shift d 



6P1/2 lifetime 6 
6P3/2 lifetime 6 



a 

P s 
6S HFS h 
7S HFS 1 
6P1/2 HFS j 
7P 1/2 HFS k 



(7P || D || 7S) 
(65 11 D || 6P 1/2 ) 
<6S||£>||6P 3/2 > 
(75 || D || 6P 1/2 ), and 
(7S\\D\\6P 3/2 ) 
same as a 
ipes(r = 0) 
Tp7s(r = 0) 

(l/r 3 ) 6 P 
(l/r 3 ) 7P 



-3.4[19] 
-4.2[-8] 
-2.6[-41] 



1.8 
-6.0 
-6.1 
-7.1 



-0.7[22] 

4.3[1] 
7.9[-31] 

-1.4 
-0.8 
-3.1 
-3.4 
2.6 
-1.5 



1.0[4] 
1.0 [43] 
2.3[22] 

3.2 
3.0 

0.2 
0.2 
0.5 



a The value for fcpNc of Dzuba, et al. is obtained using "energy rescaling" so we have used the 
corresponding "rescaled" values in the table for consistency. Blundell et al. do not rescale 
/cpnc and so we use their pure ab inito values in the table. b Refs. |fl3|JT4| ; 6 Refs. |fl2|,|6|; 
d Ref. P|; 6 Ref. [18|; f Using present work's value of (3 and a/ (3 from Ref. |L9|; g Present 
work; h Denned; 1 Ref. 0; j Ref. |2|; k Ref. g|]; 
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FIG. 1. Sample data comparing scans with and without an applied electric field. Open circles 
are with E = 707 V/cm, and the scale on the right. Closed circles are with E = V/cm and the 
scale on the left. The two lines are offset from one another and have different widths because of 
the different sensitivities to ac Stark shifts for the Ml and El transitions. 
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